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t — ' Abstract: We establish that by parameterizing the configuration space of a one- 



dimensional quantum system by polynomial invariants of q-deformed Coxeter groups it 



is possible to construct exactly solvable models of Calogero type. We adopt the previ- 



o 

dimensional representation spaces of a Lie algebra. We present explicitly the G^-case for 



ously introduced notion of solvability which consists of relating the Hamiltonian to finite 



which we construct the potentials by means of suitable gauge transformations. 



1. Introduction 



One of the ultimate goals in the study of quantum mechanical systems is to find explicit 
and possibly exact solutions for the eigensystem of Hamiltonian systems. The Calogero 
U g I and Sutherland §, |, |, § models are some of the well known examples for 
theories which are integrable and can be solved exactly, classically as well as quantum 
mechanically. The integrability of the models was established more systematically by 
relating them to Lie algebraic structures, in the so-called Hamiltonian reduction method 
H, H|, or by formulating Lax pairs and zero curvature conditions @, ||, 0, |5[ ||] . 



Relatively recent |T7], 18, 19, 20] the procedure to establish their exact solvability (which 



is conceptionally different from integrability) was put on a more systematic ground by 
relating first the coordinates of the configuration space of the Hamiltonians to invariant 
polynomials. It was shown that the differential operators in these polynomials form a 
representation for certain algebras, albeit not uniquely. Having an algebraic version of the 
model, solvability can be established thereafter by noting that the eigenfunctions form a 
flag which coincides with the finite dimensional representation space of a gl(N)-Lie algebra. 
This approach has turned out to be successful in many cases and could even be extended 



to theories which are supersymmetric [21 
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In a sequence of publications |23|, |24], 25] this procedure was reversed. Instead 
of starting with a concrete potential for a theory, Haschke and Riihl proposed to start 
with a Hamiltonian already formulated in terms of invariant polynomials and construct 
the potential from it. Hence, in this approach the solvability is already build in from the 
very start and the question is addressed if possibly one obtains new types of potentials 
which are related to solvable models, which are potentially also integrable. For several 



examples of models formulated in terms of invariants of the Weyl Group p3, E3, 24] and 



even for Coxeter groups which are not Weyl groups [25, 26 1 it was shown that this is indeed 
possible. 

The main purpose of this paper is to demonstrate that this procedure can also be 
carried out successfully for models which are related to q-deformed Coxeter groups. We 
demonstrate this for G\. At the same time we show that also for these groups the associ- 
ated Hamiltonians can be formulated in terms of the gl(N)-Lie algebra generators, hence 
guaranteeing their solvability. 

We shall focus here mainly on the construction of new potentials of Calogero type, 
adopting to a large extent the point of view of the aforementioned papers. The obvious 
question of solving the associated Schrodinger problem similar as it has been done for the 
An case in [2], |3| shall not be our concern here. While this is an interesting problem 
for future work, it appears that it is still open even for almost all of the non-deformed Lie 
algebras other than the ^Ar-series. 

Our manuscript is organized as follows: In the next section we recall the notion of 
solvability based on the fact that certain types of Hamiltonians can be formulated in terms 
of the generators of the Borel subalgebra of the gl(N)-Lie algebra. We show how from 
this formulation one may systematically construct potentials. In section 3 we assemble the 
main mathematical properties about polynomial invariants of the Coxeter group, which 
play the crucial role of coordinates in this context. In section 4 we extend these ideas to 
the q-deformed Coxeter groups. In section 5 we discuss how certain choices of the pre- 
potential lead to Cologero type potentials. In sections 6 and 7 we discuss the Calogero 
model for G2 and its q-deformed version, respectively, deriving some explicit Calogero type 
potentials. We state our conclusions in section 8. 

2. Construction of exactly solvable potentials 

We start by recalling the notion of exact solvability as proposed originally by A. Turbiner 
[17] about ten years ago. For this we require polynomial spaces of the form 



K = span<J/ 2 fc2 / 3 fc3 .../^ 



The I{ constitute some generic set of variables not further specified at this point. Evidently, 
these spaces are embedded into each other Vq C V\ C V2 C . . ., hence forming an infinite 
flag. A Hamiltonian operator TC acting on such spaces and respecting 

H: V n >-> V n (2.2) 
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possesses an infinite family of polynomial eigenfunctions. Therefore, it is natural to refer 
to such type of Hamiltonians as exactly solvable. 

It is now a matter of identifying the spaces V n , which of course allows for numerous 
solutions. It was noticed that many known models can be fitted into this scheme when 
one identifies V n with a finite dimensional representation space of a gl(N)-Lie algebra. A 
simple representation of this algebra in terms of first order differential operators is found 
when expressing the usual gl(N)-generators Eij as 



N 



= J?- = Iidj, E i0 = J i = dj, E 0i = = K Ij - Iihdk for k G 



(2.3) 



fc=2 



where k is an arbitrary constant and di = d/dli. It is easy to check that these differential 
operators satisfy indeed the usual gl(N)-commutation relations 



[Eij, Eki] = 5jkEu — 5uEkj 



(2.4) 



According to the representation ( |2.3D all Hamiltonians which are expressible in terms of 
the Borel subalgebra of gl(N), i.e. involving only the generators J° and J~, will respect 
Q2.2p and are therefore exactly solvable. Remarkably, it has turned out that many known 
solvable models of Calogero and Sutherland type can be brought into the general form 



Ti — c ij Jij Jkl + &ij J?j Jk + ^ ~^ 



(2.5) 



with cfj, c5L- 3 di, dij G M. being some coupling constants. 

Unfortunately not all models, in particular the ones we shall discuss below, can be fitted 
into the gl(N)-framework. Nonetheless, following the same ideology as outlined above, one 
can appeal to some other algebras which can be realized with different types of differential 
operators than those provided in (|2.3j). For our purposes the semi-direct sum g^M) ix ' 



will be rather useful [17|. It may be realised by the £ + 5 generators 



J 



di, 



J 



J 4 = Ifdi + ihhd 2 - Ki- 



ll 



J 3 : 



hd2 ~ Yt 

I{d 2 



for 1 < i < £, k G 



(2.6) 
(2.7) 



The further condition k G Z + guarantees that the representation is finite dimensional. 



Expressing now Hamiltonians in terms of the Borel subalgebra of g^QR) x 
J 1 for 1 < i < £ + 5 with i ^ 4, the flag space of the form 



V n = span 



|0 < h +£k 2 < n 



, i.e. the 



(2. 



will be left invariant in the sense ( |2.2j ). 

The above observations inspired the starting point of the approach in 
which is the eigenvalue equation for the function ip(I) 



23), |§ |5| 



Dip = Etp 
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with T> being a symmetric Schrodinger operator of the form 

V = ~ S d k9u d i + rkdk ■ ( 2 - 10 ) 

Here g^ 1 denotes the inverse of the curvature free symmetric Riemannian tensor = gik, 
which, in view of ( |2 . 5| ) , is at most quadratic in the coordinates I{. The functions r k are 
assumed, again in view of ( |2.5p , to be linear in the coordinates ij. In many cases these 
coordinates are taken to be invariant polynomials (for more details see below) , albeit some- 
times re-parameterizations are needed to guarantee the quadratic and linear dependence 
of gki and r&, respectively. 

Clearly, the operator ( 2.1 Q| ) is not of the usual form of a Hamiltonian, that is Laplacian 
plus potential. In order to extract a potential from this Hamiltonian one has to carry out 
a gauge transformation <p = e x ip to bring the equation (p.9[) into the more standard form 



(-A + V)V> = Eip (2.11) 
involving the Laplace-Beltrami operator in general Riemannian space 

A = iVS^ with G~ 1 = detg- 1 (2.12) 

VG k,i 

and a potential V. Extracting then from the equality e~ x T>e x = —A + V the terms of first 
and zeroth order in di, one finds 

rfc = I> i ^(2x-lnv / G) (2.13) 
V = r Ax ~ |] [dkig^dix) + QudkXdix] , (2-14) 

respectively. Multiplying now ( |2.13| ) with g^} and differentiating thereafter with d m one 
realizes that the right hand side is symmetric under the exchange m *-* I. Therefore, one 
deduces immediately for the left hand side the same symmetry 

®m ^2(9kin) = d k ^2{g m in) • (2.15) 
i i 

This equation constraints the values of and can be solved by 

rfc = (2-16) 

l 

The function p introduced at this point is referred to as pre-potential. It should be stressed 
that there is no compelling argument in this approach, which fixes this pre-potential and 
it remains subject to a convenient ansatz. Substituting ( 2.16|) back into ( |2.13| ) and ( 2.1 1 ) 
one then finds 



X 



hp + lnVG) (2.17) 



2 



V = lEA^- lE^( ln V^(mVG) -Y^UQudix) ■ (2-18) 



1 ' " 4 
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It will turn out below that the term in V which involves % 1S zero or constant. Provided that 
diQny/G) is linear in /, this would follow directly as a consequence of the assumption 

already made on r&, namely that it is linear in the variables /. In that case we can deduce 

that this term in the potential would be constant and omitting it just amounts to a constant 

shift of the ground state energy. 

Before we can specify in more detail the ansatz for the pre-potential p proposed in p4| . 

we have to gather various facts about invariant polynomials. This will make the suggested 

ansatz look very natural, albeit not entirely compelling. 



3. Polynomial invariants of the Coxeter group 

We specify here in more detail the nature of the variables / and assemble some of their 
mathematical properties. First we recall the well known fact, that to each simple root ai 
in a root system A one can associate a reflection on the hyperplane through the origin 
orthogonal to ai 

ai {£) = x-2^-^oti for 1 < i < £, x G B*. (3.1) 

These reflections constitute the Coxeter group W of rank £ or more specifically when 
2a ■ (3/(3 2 G Z for all a, f3 G A a Weyl group. One may then express each vector x G M. e as 
x = Yld=i x i a i an d associate to it a polynomial P(x±, . . . ,xg). The action of the Coxeter 
group on these polynomials is defined as 

a i P(x 1 ,...,x e ) = P(ar\x 1 ),...,a- 1 (x e )) . (3.2) 



From the defining relations ( p.l| ) and (3^2) it follows directly that by taking the simple 
roots as a basis for the action of the simple Weyl reflections acquires a particularly 
simple form 

o- i P(x 1 , ...,x e ) = P(xi, ...,Xi-i,Xi - ^2xjKji,x i+ i, ...,x e ) . (3.3) 



Here K denotes the Cartan matrix = 2cti • ctj/a'j. The special set of polynomials which 
does not change under the action of W, i.e. for which 

ail s (xi, . . . ,x e ) = I s (xi, . . . ,x £ ) for all ai G W (3.4) 

are the polynomial invariants of the Coxeter group. It turns out that a basic set of linear 
independent polynomials {I\ +Sl , . . . , Ii +Se } can be graded by the £ exponents Sj of the 
Coxeter group, with 1 < i < £. The subscripts 1 + Si indicate here the degrees of the 
polynomials. It is this set of basic invariants which one takes as the coordinates of the 
previously described Hamiltonian system. 

Let us now establish and recall some of their main properties, which we shall exploit 
below: 
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3.1 Eigenbasis of the Coxeter element 

It is clear that the choice of the basis for the coordinates will alter the form of the potential 
and a priori there is no coordinate system which is more special than another one. However, 
certain choices make the final expressions very simple and one can take here the search 
for simplicity as a guiding principle. A particularly suitable choice is the eigenbasis of the 
Coxeter element. We will see that in this basis the expressions for the polynomial invariants 
simplify considerably. 

Adopting the notations of p7| (see also references therein), we first define the Coxeter 
element a in terms of the two special elements of the Weyl group 

a± := n °* ' ( 3 - 5 ) 

as a := c_cr+. Here we have partitioned the set of simple roots into two disjoint sets, say 
a k £ <3?+ and (3 k £ <&_, by associating the values q = ±1 to the vertices i of the Dynkin 
diagram of the Lie algebra, in such a way that no two vertices related to the same set are 
linked together. The eigensystem of the Coxeter element can then be brought into the form 

aVj = e*sr'ivj and Vj = e~^ £ jk a k + ^ £ jk (3 k , (3.6) 

k k 

where we denote by £ the matrix of left eigenvectors of the Cartan matrix, i.e. 

t 

^^• fc = 4sin 2 ^ ifc , (3.7) 

3=1 

and the s, are the aforementioned exponents. Then we implicitly define a variable substi- 
tution {xi} — ► {wi} by the basis transformation 

s= y~] Xioti = ^2 dj w j a i = ^2 wiVi ■ ( 3 - 8 ) 

i i,j i 

with 

_ / e~^%- fe , for a k G 

Qk3 I i ]k , for P k e *_ ■ (3 ' 9) 



Defining then also polynomials in these new variables, we obtain as a consequence of ( |3.6j ) 
the action of the Coxeter element on these polynomials 

aP(wi, Wfj = P{w\e h , h 2 } ...,w£e h). (3.10) 

Recall here that s« + s^_j = h for 1 < i < t. Since the Coxeter element is build from 
simple Weyl reflections it follows from ( |3.4| ) that the invariants of the Coxeter group are 
also invariant under the action of the Coxeter element 



al s (x 1 ,...,x i ) = l s (x 1 ,...,x e ) =Is[^2 Cii^i, .-.,>; CuWi | • (3-11) 



To be able to compute the action of the Weyl reflections we have to express the polynomials 
in terms of the x-variables. Nonetheless, by means of (|3.8| ) we can also translate the action 
of the Weyl reflections in the x-variables to an action in terms of the w;-variables which 
allows for a more concise and possibly generic formulation of the invariants (see below). 



- 6 - 



Exactly solvable potentials for q-deformed Coxeter groups 



3.2 Universal formulae for invariants 

It is a natural question to ask, whether there exist formulae which express the invariants 
in a universal fashion, that is valid for all algebras. Indeed for invariants of degree 2 this 
is possible and we find 

I 

h{x) = yjjXj ■ .r,I\ ;j l j.i-j. Kijtj = KjiU, ti > 1 . (3.12) 

i=l i<j 

Here the ti denote the symmetrizers of the Cartan matrix, which could be avoided in the 
above expression by absorbing them into the roots, which amounts to taking the simple 
co-roots instead of the simple roots as a basis. For higher degrees we did not succeed to 
find universal formulae for the invariants. 

Changing, however, to the eigenbasis of the Coxeter element it is far more obvious how 



to write down a universal expression. From ( 3.10 ) it is clear that any invariant has to be 
of the form 

s 

I s (w) = c s (ai, . . . ,ae.)w^ . . .w^ e (3.13) 

ai ,...,a,£=l 

where the constants c s (ai, . . . , ai) are constrained as 

c s (a 1 ,...,a e ) = {^° if ^ flt = '' Etx^i = nh, n G Z _ 
I = otherwise 

Consequently, this means for instance that the quadratic invariant has to be of the form 

h = J y2c2(a i ,a£_i +1 )wiW£_ i+1 . (3.15) 

i 

One can proceed similarly for higher degrees, but it is then less obvious how to fix the 
constants c s (ai, . . . , ai). Hence, for the time being we have to rely on case-by-case studies, 
but even for explicit algebras the most generic expressions are difficult to find in the 
literature. See p8| for a complete list. 



3.3 Jacobians and Riemannians 

It will turn out that a key quantity in this scheme is the Jacobian determinant related to 
the polynomial invariants 

dl 

J = det (j) with j kl = ^±±^ . (3.16) 

dui 

The determinant J is known to possess various important properties, see e.g. 2£]: 



i) The polynomials Ii+ S1 , • • • , Ii+s t in ui, . . . ,v,£ are algebraically independent if and only 
if J / 0. 
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ii) Defining for each root a a linear polynomial 



Pa (til, ...,Ut) =J2la )u i 
i=l 



with 7$ G 



(3.17) 



such that p a (ui, ■ ■ ■ , tt^) = defines the hyperplane through the origin orthogonal to 
a. Then one can factorize J as 



J = |J |] p a (ui,... ,U£) with fl £ 



(3.18) 



where A + denotes the set of positive roots, 
iii) Defining the inverse Riemannian in terms of the basic invariants {ii+ Sl , . . . , Ii+s e } 

dhdh 



9kl = Yj 



i=l 



du; dui ' 



(3.19) 



the Jacobian determinant is related to the determinant of the inverse Riemannian as 

J 2 = G~ l = deter 1 . (3.20) 



The factorization properties ( 3.18j ) and ( 3.20 ) for J and G , respectively, will make the 
ansatz for the pre-potential p appear very natural. 

It is worth pointing out that the choice of the metric ( 3.19j ) guarantees that the Lapla- 
cian in the variables u is flat. This is easily seen by considering the change of the Euclidean 



metric tensor, i.e. g(u) 
to multiply 



under a coordinate transformation. For this we just have 



>, du m du 11 \ 
9^ = 2Z^r^f-9{u) m n = Y, 



du m du m 
dli dlj ' 



(3.21) 



with ( p. 19 ). This choice avoids the entire analysis which is needed in this approach to 



guarantee the flatness of the Laplacian as carried out in [22| 



4. Polynomial invariants of the q-deformed Coxeter group 

We extend now the previous discussion and seek polynomials which are invariant under 
(/-deformed Weyl reflections. We adopt here the notation of p^| , (see also p0| ), for more 
details on g-deformed Weyl reflections and the general context in which they emerged. 
When acting on a simple root they are defined as 

crf(aj) = aj - (25ij - [Iji] q )ai (4.1) 

with I = 2 — K being the incidence matrix of some Lie algebra and [n] q = {q n — q~ n ) /(q — 
q^ 1 ) being the standard notation for a g-deformed integer. According to the notions out- 
lined in the previous section the invariants are obviously defined by 

ajl s (xi, . . . ,X£) = I s (xi, ■ ■ ■ ,xt) for 1 < i < £ . (4.2) 
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Once again for polynomials of degree 2 we can write down a universal formula 

I 

4 = I>W +Y,n: i [K ij ] q [t j ] q x j (4.3) 

i=l i<j 

which is invariant under the g-deformed Weyl reflections. As in the non-deformed case 
higher invariants are difficult to write down in a universal form. Note also that the q- 
deformed invariants in the sense of Q4.2j) are not invariants of the (/-deformed Coxeter 



elements as defined in [3C, |27f l, as the latter does not just consist of transformations af. 



However, we can alter this definition slightly to achieve this, see section 7. 



5. Ansatz for the pre-potential 

Having fixed a set of basic invariants . . . , Ii+ Se } one assumes that the wavefunctions 

in fl2.9Q and ( 2.11[ ) depend on these coordinates, that is (p — > ip(I) and ifi — > Naturally 



one can then also view the potential as a function of the invariants, i.e. V — > V(I) and 
understand that d k = d/dlk- Defining the inverse Riemannian as in ( gig) and using 
( 3.20 ), we can re- write the potential ( |2 . 1 8| ) in the form 



V = \ Yl dudkpdip -JJ2YI 9 k id k JdiJ - i d * 9h ( d iP ~ J^ J ) 



(5.1) 



For the above mentioned reason the last term will usually drop out. It is apparent from 
this formulation, that re-scaling J by a constant will not alter the potential, a fact which 
is important with regard to the occurrence of possible coupling constants. To be more 



specific about the potential one has to choose a suitable pre-potential. In |23|, 24, |25[ the 
following ansatz was proposed 



J2lilnPi(l). (5.2) 



i=0 



The P%(I) for 1 < i < I are defined by the factorization of the determinant of the inverse 
Riemannian 

1 

J 2 = G -i = detg -i = Y[p i (r ) . (5.3) 

1=1 

Evidently, this ansatz ( |5.3| ) is inspired by the properties ii) and iii). However, there is an 
additional significant constraint namely that the Pi{I) are functions of the invariants /, 
which as was argued above is needed to guarantee the solvability. Often one would like to 
obtain also an additional harmonic confining term proportional to uf in V . This is 

easily achieved by including also a factor of the form Pq ~ exp(^^ uf) into the ansatz for 
the pre-potential (|5.2[) . For the reasons outlined in section 2, the entire Hamiltonian, that 
means also this term, has to be expressed in terms of invariant polynomials. Usually we 
can take Pq = exp^)- 
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Substituting the ansatz ( |5.3[) into ( |2,18| ), the potential acquires the form 

V = I E U - i) ^(mP^On^) - ~ £ (7* - 5) & |>H ftOnPi)] • (5-4) 
*• - fe ' ^ ' i,k,l ^ 7 



i,j,k,l 



As was pointed out in p4|], it will turn out that the terms with i 7^ j in the first term 
are constants (even zero) most of the time, which therefore can be dropped safely by just 
shifting the ground state. The other motivation for the ansatz fl5.3| ) is that one would like 
the two terms in ( |5.l|) or ( [2.18D to combine naturally. 

Comparing now ( PP|) , (gl8j) and (g!20|) we can proceed and exploit the fact that the 
Pi factorize further into linear polynomials 



Pi = n t^)] 2 • 



(5.5) 



aSA 



The relation (|5.5|) is here the defining relation for the set of positive roots A$ . Exploit 



ing the fact that the p a (u) are linear in u, see ( |3 . 1 7| ) , and changing from the invariant 
polynomials as coordinates to the u- variables, we obtain 



V 



E 



1 -j 



7<7i 



E 



(fc).(fc) 



aeA^,/3eA 



2^fc7o 7/3 

Pa(u)Pf3(u) 



(5.6) 



Recall that the 7^ are defined in (3.17). The expression Q5.6| ) constitutes a general for- 
mula for potentials when starting with any Coxeter group. This structure will survive 
in different coordinate systems, as clearly the p a (u) remain linear after linear coordinate 
transformations. Often one changes the coordinate system by using an explicit represen- 
tations for the roots in an orthogonal basis, which one may find in various places of the 
literature, e.g. I^Sj] . Then (|5.6|) enables one to write down directly the potentials associated 
to any Coxeter group. In practice it turns out that only the diagonal terms in the sum 
survive, i.e. a = (3, such that the potentials will always be of Calogero type. We omitted 
in ( |5.6| ) the term resulting from the last term in (5.4) and also a possible constant. 



6. Prom G2-polynomial invariants to the Calogero models 



In [19] it was shown that the Calogero models for three particles exhibit an underlying 
(^-structure, which can be exploited to establish their solvability. In [24] this procedure 
was reversed and it was shown that the approach outlined in section 2 and 3 indeed yields 
potentials of the Calogero type when one starts with a G2-structure. Here we recall briefly 
the procedure, mainly to set the scene for the q-deformed treatment below, but also to 



establish a few facts not pointed out so far. In particular, in [19| as well as in [24| not the 
most general GVinvariants were used. As we will show below the most generic invariants 
involve some arbitrary constants. The obvious question to ask is whether one obtains a 
new type of potential when using the procedure outlined above in terms of these generic 
coordinates, possibly involving additional coupling constants. 
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Let us start with the action of the GVWeyl reflections on the simple roots 

c"i(ai) = —a\ o~i{(*2) = 3ai + c*2 (6-1) 

0"2(«i) = «i + «2 cr 2 (a2) = — «2 • (6.2) 

Then for a general vector x = x\a\ + X2OL2 m we have 

o~i(x) = (3x2 - xi)ai + x 2 «2 (6.3) 

a 2 (x) = x\a\ + (x\ - x 2 )a 2 , (6.4) 

such that 

(7iP(xi,x 2 ) = -P(3x 2 - xi,x 2 ) (6.5) 

£72-P(xi,X 2 ) = P(xi,Xi-X 2 ) • (6.6) 



The equations (|6.5[) and (|6.6[) follow also directly from (|3T3| 



Using ( |6.5| ) and ( |6.6[) we can now generate solutions to the equation (3^), i.e. construct 
the invariant polynomials. The procedure is straightforward. We simply write down the 
most generic expression for a potential candidate for a polynomial invariant I s (x) of degree 
s similar to the form as in ( 3.13j ) with arbitrary coefficients c s , but now in terms of the x- 



variables. Acting then with all simple Weyl reflections on this polynomial and demanding 
invariance (3.4) leads to a system of equations which determine the c s . Depending on the 



degree and the algebra this might not yield enough equations to fix all constants and one 
ends up with expressions still involving free parameters. In this manner we find as generic 
invariants 

h = k 2 Qxf + x\ - xix 2 ^ (6.7) 

-—x? + -xfx 2 - -x\x\ + hx\x\ - §x x x\ + 2xl)+ ^2 ■ ( 6 - 8 ) 

Here K2, kq and kq remain arbitrary constants. We see that besides an overall constant, 
which is naturally always present, Iq also involves an additional free parameter k§. 

A first restriction on possible values the constants might take comes from the fact that 
we want I2 and 1% to be algebraically independent. To establish this we compute first the 
Jacobian determinant for these invariants according to the definition (|3.16j ) 

2 

J = -K2KqX\X2{xi - x 2 )(xi - 2x 2 )(xi - 3x 2 )(2xi - 3x 2 ) . (6.9) 

As we expect from property i) stated in section 3 and the explicit expression in (|T^) we have 
to keep kq 7^ in order to guarantee the algebraic independence of I2 and Iq. Obviously, 
for kq = we have Iq = HqI^. We also note that kq can remain completely arbitrary in 
this context, but as we see below, we can not simply set it to zero for our purposes. 

Alternatively, we can compute the Jacobian determinant by an entirely different for- 
mula, namely ( p,18| ), and thus confirm the computation which led to ( |6.9| ). To be able to 



use ( 3.18 ), we recall that the positive roots of G2 are, (see e.g. [31, p9j])) 

A, 2 = {ai, a 2 , ai + a 2 , 2«i + a 2 , 3ai + a 2 , 3ai + 2a 2 } , (6.10) 
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such that the hyperplanes result to 

p ai = 2xi ~ 3x 2 , Pa x +a 2 =X\ - 3x 2 , P2 ai +a 2 = x l> /g 1; q 

P3ai+Q 2 = x l - x 2, Pa 2 =X\-2x 2 , P3 ai +2a 2 = X\. 

Assembling the p a into the product ( |3.18|) the result confirms our findings ( |6.9|) with 
[i = (2/3)k 2 kq. Notice, that the constants k 2 and k% organize in a separate factor in J, 
such that different choices, apart from k 2 = or kq = 0, will not alter the polynomial 
structure of J. From ( pM| ) we deduce that in the second term of the potential the overall 
factor in J just cancels, such that the constants K{ completely drop out from this term. 

So far, we have only two coordinates. To incorporate a three body interaction we need 
one more coordinate. Let us therefore choose an orthogonal basis in M 3 for the simple 
roots «i = E\— e 2 and a 2 = — 2ei+ £2+ £3, with e% • Ej = dij (see e.g. Then we can 

introduce a new set of variables via the relation 



x 



= xiai + x 2 a 2 = (xi - 2x 2 )ei + (x 2 - x\)e 2 + x 2 e 3 = y\e\ + y 2 e 2 + y 3 e 3 , (6.12) 



with the built-in constraint y\ + y 2 + 7/3 = 0. In these variables the action of the Weyl 
group on the same polynomials becomes 

0-1^(2/1,2/2,2/3) = ^(2/2,2/1,2/3) (6.13) 
o 2 P(2/i,2/2,2/3) = P(~yi, -2/3, -2/2) • (6.14) 

In principle, we could proceed as above, i.e. writing down generic expression with arbi- 
trary coefficients and use directly ( |6.13| ) and ( 6.14 ) to find invariant polynomials in the 



y-variables. However, since we have increased the number of coordinates also the amount 
of unknown coefficients grows and we will end up with polynomials involving many more 
free constants than just the three we are left with when using the x-variables. Instead, as 



we know the invariants already, we can also use in ( |6.7| ) and (6.8) directly the substitutions 
x\ — > —y\ — 2y 2 and x 2 — > —y\ — y 2 , such that we obtain 

h = K 2 (yj + y\ + 2/i2/2)/3 (6.15) 
h = 27 («6 - 2K6)(yf + + g(«6 - 2K 6 )(yi2/2 + 2/12/2) 

+ i(5 K6 + 2^ 6 )(yfy 2 2 + yfyf) + ^(40 K6 + 7K 6 )y 3 y 3 . ( 6 .16) 



Only for the special choice k 2 = — 3, kq = 2kq and kq = 1 the expressions (|6~15|) and (gig ) 



reduce to the coordinates Ai, A2 used in [19]. Notice, that in the y-variables the invariants 



become symmetric polynomials [32] 



In order to reproduce the Calogero potentials we need to make yet another coordinate 
transformation and introduce Jacobi relative coordinates 

1 3 

2/i = Zi~ -^_j z i » ( 6 - 17 ) 
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which separate off the center of mass motion. The constraint in the y-coordinates is now 
replaced by z\ + zi + 23 = 3Z, where Z is constant. In these coordinates we compute the 



inverse Riemannian (3.19) to 



3 



9ij 



S ^dI 1 dI l = K / 2/3/ 2 2/ 6 
^ «9z fc <9z fc ^ 2/ 6 4 [2^6/6 + / 2 3 (4^i - Rl)/4] ) . : 



Apparently <^ is not quadratic in the variables /, which as we discussed is a necessary 
requirement to be able to bring the Hamiltonian into the form fl2.5| ) and hence ensuring 
the solvability of the model. However, for R,q = 2kq one can choose a different set of 
variables 1% = t%, Iq = r|, see [19|, such that dig^ dj is of the desired form in the 



r variables. Alternatively, one can also use the representation of a different algebra to 
establish solvability Let us from now on take kq = 2kq and K2 = 1 but leaving kq 

arbitrary. Then we compute from ( |6.18| ) 



Cr 1 = detrT i = 4/ 6 (4k 6 /| - / 6 ) . (6.19) 



For the pre-potential we make now an ansatz according to ( |5.3j ), where we also include the 
previously mentioned Po-term 

P = e h , Pi = 4k 6 / 2 3 - / 6 , P 2 = / 6 . (6.20) 

From this we compute with formula (^J) the potential to 

V = -§-h + 3A lK6 — -3— - + 3X 2 kS (6.21) 
6 4k 6 /| - / 6 h 

fc=l 1<«<J<3 J 1<*<J<3 J 

where the coupling constants are to = 7 /(3y / 2) and A, = 27? — 1/2 for i = 1,2. The 
potential in the form ( |6.22 ) corresponds to so-called rational GVmodel [33], which reduces 



to the Calogero model JT], ||, |j| when the three-particle interaction is switched off, i.e. for 
A2 — ► 0. Notice that the coupling constants j , 7 X and 7 2 which enter the scheme through 
the ansatz for the pre-potential just reparameterize the coupling constants of the GVmodel. 
Note also that the constant kq has dropped out completely, such that any choice, apart 
from kq = 0, will yield the same potential ( |6.22 ). Hence, the ambiguity in the choice of 



the invariant polynomials as coordinates has no bearing on the physics. 

The invariants acquire a particularly simple form when we use the eigenbasis of the 
Coxeter element. The transformations outlined above yield for the G2-case 

Xl = V3{e- i7r/6 Wl - e- in5/6 w 2 ), x 2 =w 1 +w 2 , (6.23) 

such that the invariants simplify considerably 

I2 = K2W1W2 and Iq = Ke(wx + luf) + Kewfuu^ ■ (6-24) 
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We then find for the Jacobian 

J = -§k 2 k§{w\ - wl) (6.25) 

= -§K2Kq{wi + W2)(W\ - W2){w{ - WiW 2 + w\)(w\ + WiW 2 + w\) , (6.26) 

which could of course be used to construct the potential in these variables. As this type of 
factorization of J involves quadratic polynomials we will end up with potentials not quite 
of Calogero type. To achieve this we would have to factorize the last two terms further 
involving complex coefficients, but in that case the individual two particle interactions 
terms would be complex. Solvability is only guaranteed when we can express the factors 
in terms of the invariant polynomials. 

7. Exactly solvable potentials from q-deformed G 2 -polynomial invariants 

Let us now extend the previous analysis to the q-deformed case. To commence we need 
to evaluate the q-deformed Weyl reflections a\ as defined in fl4.1| ) for which we require the 
q-deformed Cartan matrix. In our conventions it reads for the G2-case 

With K q at hand we can now seek invariant polynomials according to the definitions ( f4.1| ) 
and Q4.2p . We proceed in the same manner as for the non-deformed case and start with 
generic expressions I$(x) for polynomials of degree s as in ( 3.13 ) and fix the constants as 



outlined above. For generic deformation parameters q we did not find invariants. However, 
if we parameterize the </'s as 



2 1 f-i , n 2lT \ If-, , o„„ 2lT \ 2 



q = h i+2cos tJ + vv 1+2cos ^J " 4 (7 - 2) 

with h being some integer, the q-deformed Cartan matrix becomes 

Implicitly, we used here ( [7.3| ) to define some q-deformed roots af. Clearly for h = 3 we 
recover the Cartan matrix of A2, for h = 4 we obtain the one of C2 and h = 6 corresponds 
to G2. For the values (|7.2j ) of q we find there exist always the invariants l\ and I?. From 
the above mentioned arguments this suggests that the exponents of this algebra G\ are 1 
and h — 1 . This assertion is supported by the observation that the formula for eigenvalues 
of the Cartan matrix ( p.7|) still holds for the q-deformed case (|7.3j ) when taking s\ = 1 and 
s 2 = h—1. 

As we argued in the previous section, it is difficult to find generic expressions for the in- 
variants in the x-variables. However, as we will see in the eigenbasis of the Coxeter element 
this task simplifies drastically. According to (|3.8|) and (|3.9|) we have the transformations 



l + e~— l + e— \ _ i [ 1-1- e~ \ 

\w w = tt- x , (7.4) 

1 1 J 2sin^ I -1 l + e— r J 
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such that the q-deformed Weyl reflections in the iu-variables simplify to 

o-\{w\) = w 2 , o-\(w 2 ) = wi, o-\(w\) = e~w 2 , o-\(w 2 ) = e ^~wi . (7.5) 

The o"|-transformations dictate that the invariants have to be symmetric in w\, w 2 and the 
^^-transformations constrain their overall degree. With ( |7.5| ) we can easily find the most 
generic expressions for the invariants 

If = k 2Wi w 2 and I q h = Khivft + u#) + R h {w 1 w 2 ) h/2 ■ (7.6) 

where k^ = for h being an odd integer. We can now transform back to the x-variables 
and confirm for instance the generic formula ( |4.3| ) for the invariant of degree 2, which still 
takes on a fairly simple form 

I9 = Jrr — -— ( x? + 4cos 2 — Xn — 4cos 2 — ) . (7.7) 

2 4sin 2 (27r//i) V 1 h 2 h 1 

On the other hand, the expressions for the I? are already quite cumbersome, albeit it is 

clear how to construct them from ( |7.6| ) and ([7.4|). 

As we saw in the previous section it was crucial to change the coordinate system yet 

further to recover the Calogero potentials in the usual form. We proceed here similarly. 

Let us choose first an orthogonal basis for the two simple q-deformed roots in R 3 

a r~ 27T 27T 2ir 2ir <— 2tt , , /— 

a\ = (V3cos— — h sin — — , — 2 sin — , sin — V3cos— — )/v3 (7-8) 

h h h h h 

c& = (—Vs — V3 cos — sin — , 2 sin — — , y/3 + v3 cos — sin — — )/VS . (7.9) 

h h h h h 

The inner products of these roots are a\ ■ a\ = 2, ■ = 8 cos 2 ^ and a[-a\ = —4 cos 2 ¥ 

such that we recover the q-deformed Cartan matrix according to Q7.3j ). Of course the 

choices ( |7.8|) and (17.9| ) are not uniquely determined. As an additional selection criterion 

we demand that If will be of an analogous form ( |6.15| ) as in the non-deformed case for 

all choices of h, such that it will be ensured that we can express in terms of If. At 

the same time this will ensure that lb = ex P(lf ) yields the harmonic confining potential, 

similarly as for the standard G2-case. The above choice for the simple roots induces a 

definition for new variables 

y = xia\ + x 2 a q 2 , (7-10) 
which satisfy the constraint y% + y 2 + 7/3 = 0. In turn this means that we can replace in 

such that If is indeed always of the form ( |6.15| ) 

- i^km w + yl + M2) " sfej (t 4 - 3Z2 ) ■ (7 ' 12) 

Again it is obvious how to obtain the expressions for 1^, but they turn out to be more 
cumbersome. From (|7.6| ) it is apparent that the case of even and odd Coxeter number 
exhibit different behaviour. We treat them now separately and supply for each case an 
explicit example. 
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7.1 Even Coxeter numbers, h=8 



As the obvious difference between the odd and even case we found that in the even case 
the additional constant kh is entering the procedure. We shall see that there are also other 
more profound differences. We present the case h = 8 in detail. For this the relation (fj 
simply yields 



1 



-(l + \/2) + V 2 ^-! 



and 



[3], = 2 + V2. 



(7.13) 



We can now take the invariants as given by the relations ( |7.6[ ) and carry out the substitu- 
tions w — > x — > y — > z specified above. In terms of the Jacobian relative coordinates Zi the 



inverse Riemannian (3.19) results to 

h 



9i 



^2 dIi dI 3 



4/« 



k=l 



dzk dz k 



«2 



4/ 8 mi/4 [2k 8 i 8 + i$(~4-44)/4 



(7.14) 



At this point we are facing a similar problem as in the non-deformed case, that is 
h = 6, namely that is not of degree two in the variables I{. Whereas for h = 6 one may 
find a suitable variable transformation, we did not succeed in this case. Nonetheless, the 
solvability of the model may now be guaranteed now by relating the model to gl2(K) x M 4 , 
rather than gl(N), see (|J) and (fD|). 

Keeping from now on k 8 = 2n% and also K2 = 1, we can bring the Hamiltonian into 
the desired form for the Schrodinger operator Q2.5| ), ( |2.10| ) 



V 



-J 2 J 1 - 8J 3 J X - 64k 8 J 3 J 8 + 



-16^8 



;i-7a) + 



9 



4(7i +7 2 ) -5 



11 



K 8 



J 



(7.15) 



To turn this operator into the standard form ( [2.11 ) we follow the procedure outlined in 
section 2. First, we compute 



G 



-1 



det^" 1 = -16/ 8 (/ 8 



4k 8 J|) 



(7.16) 



Before we proceed further to compute the potential, let us see if we still have a relation 
between J and G _1 of the type ( p.20|) for the q-deformed algebra. In particular, we wish 
to see whether a relation of the type ( |3,18| ) still holds. For this purpose we need first of 
all a notion of positive q-deformed roots. We assume that these roots are generated in a 
similar way as the ordinary roots, i.e. by repeated action of the Coxeter element. Defining 
a q-deformed version of this 1 

(7.17) 



1^1 



2 U 1 



1 This q-deformed Coxeter element differs slightly from the one denned in |30| , EtJ, as here there is no 
T-transformation involved. The Coxeter element in [BOl |27j| is only of order h up to some factors of q. 
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we compute the entire set of q-deformed roots A q lying in the orbits of the simple q- 
deformed roots ct\, a l q (a q ), a q (a q ), . . . 





a\ 


a\ 


°\ 


-a\ = -{a\ + a q 2 ) 


a% = (2 + V2)a\ + (1 + V2)a q 


-I 


-a\ = -(1 + V2)a\ - V2a q 2 


a q 7 = (2 + 2V2)a{ + (1 + s/2)a q 2 


°\ 


-a\ = -(1 + V2)al - ol\ 


a q 8 = (2 + \ft)a\ + a\ 




—a\ 


-(A 



(7.18) 



Note that the order of the Coxeter element ( [7,17 ) is indeed h = 8, i.e. a* = 1. We adopt 
now the same notion for positive and negative roots as in the non-deformed case, that is 
we call a = n iOt q a positive root if all coefficients rii are positive. With this notion the 
set of the 2h roots can be separated equally into h positive and h negative roots. We have 
verified this statement up to h = 20 of even Coxeter numbers, which strongly suggests that 
it holds in general. Using now (7.18) we can compute the hyperplanes through the origin 
to all positive q-deformed roots 

Pa q 2 =Xi- 2x 2 

p a% = V2 Xl - 2(1 + s/2)x 2 

Xl 

(2 + V2)<ci - 2(1 + y/2)x 2 ■ 



p al = 2xi - (2 + V2)x 2 
p a% = V2~ Xl - (2 + y/2)x 2 



Pal 



X2 
Xi 



Pal =X X -X 2 



Pa q 7 
Pal 



(7.19) 



Comparing now with (fnj) we have once again a relation between J and G 1 of the type 
( 3.20 ) where J can be expressed as a product of hyperplanes ( 3.18 ) 



G- 1 



.-2 



4(3 + 2^2) 



][ (PatY 



(7.20) 



aieA2 



The two factors in ( 7.16 ) admit yet a further interpretation. Organizing the roots into two 
sets Af and A? of short and long roots, respectively, we find the identities 



-f n (p^ 2 = /s 



4/ig/| and 



Kg 



16 ( 3 + a«eA t 

According to (|5.3|) we make now the following ansatz for the pre-potential 



(7.21) 



o 



Pi 



(7.22) 



From formula (|5.4| ) and including also the Po _ t erm we then compute the potential to 



v- 7 «r 



16AlK8 



T 3 



I 3 

- + I6A2K8Y- 



(7.23) 



with X{ = (7? — 1/4) for i = 1,2. Using the above mentioned identities or directly (|5.6|), 
we can also re- write this potential in terms of the z- variables. First of all we compute 
1 „ „ ^ r i2 



Pi 



4233 
1 

4334 



(zi - z 3 ) 2 (zi + z 3 - 2z 3 ) 2 11 [(1 + eV3) Zl + (1 - eV3)z 3 - 2z 2 



e=±l 



Yl [(l-eV3-sVQ)z 1 + (l+sVs + eV6)z 3 -2z 2 



e,e=±l 



(7.24) 
(7.25) 
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We find 

Ai 3Ai 



2 " t^l' K ' (*l-*3) a ' (zi+z 3 -2z 2 ) 2 

+ E — 2 (7-26) 

«5hi [(1 + e>/3)«i + (1 - £^3>3 - 2z 2 ] 

6(2 + £gV2)A 2 

£,£±l [(1 - eV3 - e >/6)zi + (1 + e^3 + £\/6> 3 - 2z 2 ] 2 

with u> = 7 \ // 3/(2v / 2)- We omitted here a constant which contains the center of mass 
coordinate. Remarkably, all off-diagonal terms, that is terms in ( |5.4| ) with i ^ j, cancel 
each other. This potential has a very similar structure as the usual Calogero potentials 



( 6.22 ), but it involves now deformed two and three-particle interactions. We find similar 
structures for higher values of h. Responsible for this structure is the fact that we can still 
factorize J, and therefore G _1 in terms of products of hyperplanes as in ( 7.20|) . Remarkably, 
these potentials are all exactly solvable by construction. 

7.2 Odd Coxeter numbers, h=5 

The structure for theories with odd values of the Coxeter number is somewhat different. Let 
us consider h = 5 in more detail. In that case the relation for the deformation parameter 
( |7.2j ) simply yields a root of unity 

o = e^° and [3], = *±^. (7.27) 



Replacing now in equation ( |7.6| ) the variables w — > x, the invariant of degree 5 in the 
x-variables is still not too lengthy, unlike for greater values of h, and reads in this case 

if = K5 Q(3 - V5)xjx 2 - 1x\x\ + (1 + Vh)x\xl + V5)x lX {\ . (7.28) 

Obviously l\ and if are algebraically independent as one is of even and the other of odd 
degree, respectively. As for the even case, we can proceed and carry out in ( 7.28; ) the 



substitutions x — > y — > z such that the invariants are expressed in terms of the Jacobian 
relative coordinates Zj. In these coordinates the inverse Riemannian ( |3.19| ) results to 

fZi dz k 9z k \ 5if Kg//^ / .. 

Once again we have the problem that is not of degree two in the variables ij. As for 
h = 8 we can relate once more to the gl2(K) x M 5 algebra, see (|2.6| ) and (|2.7| ). From now 
on we keep k 2 = 1 and bring the Hamiltonian into the desired form (|2.5| ), (|2.10[ ) 



V = (V5-5) Q J 2 J 1 + J 1 J 3 + 54 J 5 J 9 - ( 7l - 1 - |) J 1 - I J 2 + I J 3 ) 



(7.30) 
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We proceed similarly as in the previous subsection and compute from ( [7.29 ) 

5 



G 



-i 



detg- 1 = -(V5-3)(L 



(7.31) 



In order to obtain Calogero type potentials it is vital to factorize G _1 further into linear 
polynomials. Let us proceed analogously as for even Coxeter numbers and compute the 
orbits of the q-deformed Coxeter element ( f7.17|) . We find 





q 

a\ 


q 


< 


«3 = -(<*! + «!) 


a q 7 = 1(3 + >/5)a? + 1(1 + V5)a| 




a\ = -1(1 + \/h)a\ - jCVS - l)a q 2 


ag = 1(1 + V5)al 


°\ 


a\ = |(\/5 — 1)«2 


a q 9 = -|(1 + V5)(a? + a\) 


A 


a% = 1(1 + y/l)a\ + <4 


af = — ^(3 + \/5)«i — ol\ 


A 


al 


a\ 



(7.32) 



Note that the order of the Coxeter element (7.17) is still h, that is in this case = 1. 
For odd values of the Coxeter number we can still separate the roots into positive and 
negative roots, but now the negative roots can no longer be obtained by reversing the 
signs of all positive roots. There are now positive roots without a negative counterpart. 
Unfortunately, as a consequence of this the factorization property ( 3.18| ) does no longer 
hold in its stated form. Nonetheless, one can still select some hyperplanes obtained from 



the root system (7.32) and factorize G , albeit now the selection principle does no longer 



favour the positive roots and is less clear. We compute 



Pa\ 


= xi - 1(3 + \/5)x 2 


Pal = 


Xl - 


2x 2 


Pal 


= xi - (2 + V5)x 2 


Pa q 7 = 


Xl + 


1(1 + y/b)x 2 


P^ 


= 5(1 + y/5)xi - x 2 


Pal = 


Xl - 


1(3 + V5)x 2 


Pal 


= Xl — 2X2 


Pa% = 


Xl + 


1(1 - V5)x 2 


Pal 


= - l)xi + X2 


P< = 


= xi - 


- 1(^/5 -l)x 2 



(7.33) 



and construct from this 



J 



Pa\PalPalPalPa\ ■ ( 7 - 34 ) 

Now, unlike as in the even case, the splitting into long and short roots no longer corresponds 

to factors in terms of hyperplanes. 

For the pre-potential we make now the ansatz 



Pn 



J2 



P^ = It 



(7.35) 



With formula ( ^.4[ ) we then compute the potential in terms of invariant polynomials to 

1 - ±= \h + 5(^5 - 5)A«| r2 J L r5 , (7.36) 



with A = (7! — 1/4). Once again we can also re- write V in terms of the z- variables. First 
we factorize 

2 



1 



3445 



(^i-^) 2 n 



e,e=±l 



:i + W3+-^|)z2 +g + (l 




2z 2 



(7.37) 
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from which we deduce the potential to 



1 3 A 



A ^ 6(l + e/\/5)A 



2 fc=l (Zl " Z3) e,e=±l [(1 + ^3 + ^)z 2+g + (1 - ^3 + ^)z 2 _ g - 2z 2 

(7.38) 

with cj 2 = 3(3 - \/5)7^/20. Once again all off-diagonal terms cancel each other and as in 
the even case this potential is also of Calogero type. We find similar types of potentials for 
higher values of the Coxeter number. 

8. Conclusions 

It has been shown previously that solvability of certain types of Hamiltonians can be es- 



tablished 1 17, [18|, [L9|, 20] by relating the differential operators inside the Hamiltonians, 
that is essentially the Laplace operator, to a representation of the gl(N)-Lie algebra. This 
formulation can be made very systematic by associating the differential structure to poly- 



nomial invariants of Coxeter groups. This led the authors of |2^, 23, 24, ^] to propose 
a procedure which allows to construct solvable Hamiltonians by taking the structure of 
the polynomial invariants as a starting point. Here we showed that this procedure can 
be extended successfully to polynomial invariants of q-deformed Coxeter groups. We con- 
structed some potentials resulting from these type of invariants. Due to the fact that the 
Jacobian determinant can still be factorized in terms of linear polynomials the resulting 
potentials are of Calogero type. 

There are several open issues and unanswered questions which deserve further inves- 
tigations. Clearly it would be interesting to carry out the outlined procedure explicitly 
for other algebras than G\. The presented example indicates that one can expect similar 
structures beyond G\. Eventually one should aim at a unified formulation, as opposed to 
case-by-case studies, analogously to the non-deformed CclSG clS indicated in section 5. Cru- 
cial will be here the factorization of the Jacobian determinant J. In the presented example 
this works nicely for even Coxeter numbers, but for odd h the example hints that one 
possibly has to employ a different q-deformed Coxeter transformation in order to obtain a 
definite criterion for the selection of the hyperplanes which yields the factorization of J. 

To achieve a unified formulation it will be important to have systematic and generic 



expressions for the polynomial invariants 28]. In the above analysis we have seen that 



the choice of a suitable basis is absolutely crucial for this task. The favoured one is the 
eigenbasis of the Coxeter element as we have demonstrated. 

Since we have shown that one can extend the approach from Coxeter to q-deformed 
Coxeter groups, it is also natural to suspect that one might as well employ it for reflection 
groups the general type introduced in p4| . 

Naturally it appears also possible to construct potential of Sutherland type by using 
different types of coordinates [pgj ]. 

To find the explicit wavefunction for the above Hamiltonians is now also an obvious 
question to ask. Following the quoted literature there is a straightforward procedure to 
construct them from the above mentioned results. In this context one might also address 
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technical questions as for instance self-adjointness similar as it has been done in Calogero's 
original work for hardcore boundary conditons M, 0, H], For slightly more general boundary 
conditions, see e.g. |35|| . Conceptionally one should stress that solvability in the sense 
provided here is far more constructive with regard to this question than integrability. The 
latter usually just guarantees the existence of exact solutions, whereas solvability is already 
tied closely to the explicit solutions. 
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